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ABSTRACT 

These are the proceedings of the workshop on "Lattice QCD, Chiral Perturbation 
Theory and Hadron Phenomenology" held at the European Centre for Theoretical 
Studies in Nuclear Physics and Related Areas from October 2 to 6, 2006. The 
workshop concentrated on bringing together researchers working in lattice QCD 
and chiral perturbation theory with the aim of improving our understanding how 
hadron properties can be calculated and analyzed from first principles. Included 
are a short contribution per talk. 



* This workshop was funded by the ECT* and by the EU Integrated Infrastructure Initiative 
Hadron Physics Project under contract number RII3-CT-2004-506078. Work supported in 
part by DFG (SFB/TR 16 "Subnuclear Structure of Matter" and FOR 465 "Forschergruppe 
Gitter-Hadronen-Phanomenologie" ) and by BMBF (research grant 06BN411). 



1 Introduction 



As a result of developments in lattice field theory and computer technology, the 
first full QCD calculations in the chiral regime are becoming available now. An 
essential tool for extracting hadronic quantities from lattice QCD is chiral per- 
turbation theory. Chiral perturbation theory is the effective field theory of QCD 
and allows to analyze the quark mass dependence of observables in a model- 
independent manner. Since the technologies for large scale lattice calculations 
and for detailed chiral perturbation calculations are each so different and so de- 
manding, the same theorists generally do not pursue both aspects. Hence, we 
considered it to be valuable to bring together practioners of both aspects of lat- 
tice calculations, and deepen each others' appreciation of the issues involved, thus 
opening a new era for confronting experiment with solutions of QCD from first 
principles. 

Thus, such a meeting was organized at the ECT* (Trento) from October 2-6, 
2006, with financial support from the ECT* as well as the I3HP networks N2 and 
N5 and by EURONS. The meeting had 36 participants whose names, institutes 
and email addresses are listed below. 35 of them presented results in presenta- 
tions of various lengths. Ample time was left for discussions on the various talks 
each day. In addition, there was also a special discussion session "Lattice meets 
chiral perturbation theory" convened by Hartmut Wittig, in which the interplay 
between the various issues as seen from the two communities were addressed. 
A short description of the contents of each talk and a list of the most relevant 
references can be found below. We felt that this was more appropriate a frame- 
work than full-fledged proceedings. Most results are or will soon be published 
and available on the archives, so this way we can achieve speedy publication and 
avoid duplication of results in the archives. 

Below follows first the program, then the list of participants followed by the 
abstracts of the talks. The summary of the aforementioned discussion session is 
given at the end of these proceedings. Most talks can also be obtained from the 
workshop website 

http : //www . ect . it/ 

We would like to thank the ECT*, in particular Serena Degli Avancini, for 
the excellent organization of the workshop and all participants for their valuable 
contributions. We believe that this was only the first workshop of this kind and 
look forward to similar meetings in the future. 



Ulf-G. Meifiner and Gerrit Schierholz 
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2 Program 



Monday, October 2nd 2006 

Morning session, chair: Evgeny Epelbaum 

9:45 Ulf-G. MeiBner / Opening Remarks 
Gerrit Schierholz 

9:50 Luigi Scorzato (Trento) Dynamical simulations with twisted mass 

10:25 Gerrit Schierholz Lattice QCD simulations at small quark masses 

(Hamburg) 

11:00 Coffee 

11:30 Christian Lang (Graz) Hadrons for chirally improved fermions 

12:15 Meifeng Lin (New York) Dynamical domain wall fermion simulations 

and chiral perturbation theory 

13:00 Lunch 
Afternoon session, chair: Bugra Borasoy 

14:30 Christopher Aubin Applications of staggered ChPT 
(New York) 

15:15 Stephan Diirr (Bern) t Are there alternatives to rooted staggered fermions? 

15:50 Coffee 

15:15 Sinya Aoki (Tsukuba) ChPT and lattice QCD with Wilson-type quarks 

17:05 Discussions 

18:00 End of Session 



Tuesday, October 3rd, 2006 



Morning Session, chair: Meinulf Gockeler 



09:00 Oliver Bar (Berlin) 

09:35 Ulf-G. MeiBner (Bonn) 
10:10 Veronique Bernard 
(Strasbourg) 

10:45 

11:15 Matthias Schindler (Mainz) 
11:50 Evgeny Epelbaum (Jiilich) 
12:25 

Afternoon session, chair: Christopher Aubin 



Lattice QCD with mixed actions: 
Overlap fermions on a twisted mass sea 
Thoughts on chiral extrapolations for excited states 
Chiral extrapolations of baryon properties 

Coffee 

Some topics in manifestly Lorentz-invariant CHPT 
Chiral extrapolations in few-nucleon systems 
Lunch 



14:30 Bugra Borasoy (Bonn) 
15:15 Jose A. Oiler (Murcia) 
15:50 

16:20 Hermann Krebs (Jiilich) 
16:55 Dina Alexandrou (Nicosia) 
17:40 
18:00 



Nuclear lattice simulations 

Chiral non-perturbative study of pseudoscalar self-energies 

Coffee 

Scalar two-loop diagrams on the lattice 
iV and iV to A form factors in lattice QCD 

Discussions 
End of Session 
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Wednesday, October 4th, 2006 



Morning Session, chair: Akaki Rusetsky 

09:00 David Richards Hadron structure using DWF quarks on an asqtad sea 

(Newport News) 

09:45 Wolfram Schroers (Zeuthen) Distribution amplitudes from the lattice 
10:20 Dirk Brommel (Hamburg) Pion structure from lattice QCD 

Coffee 

Nucleon structure functions and form factors 
Bk from dynamical domain wall fermion simulations 
Lunch 



10:55 

11:55 Dirk Pleiter (Zeuthen) 
12:10 Enno Scholz (Upton) 
12:45 



Afternoon Session, chair: Meifeng Lin 

14:30 Hartmut Wittig (Mainz) Low-energy constants for AS* 



15:05 Akaki Rusetsky (Bonn) 
15:40 

16:10 Meinulf Gockeler 

(Regensburg) 
16:55 Tobias Gail (Garching) 
17:30 
18:00 



1 transitions 
The A-resonance in a finite volume 

Coffee 

Lattice QCD on smaller and larger volumes 

Utilizing covariant BChPT for chiral extrapolations 

Discussions 
End of Session 



Thursday, October 5th, 2006 

Morning Session, chair: Oliver Bar 

9:00 Christoph Haefeli (Bern) Aspects of CHPT at large m s 

9:45 Poul Damgaard A chiral two-matrix theory and the chiral Lagrangian 

(Copenhagen) 
10:30 Coffee 

11:00 Jack Laiho (Batavia) The B — > D*£u form factor from lattice QCD 

11:45 Michele Delia Morte Heavy quark effective theory on the lattice: 

(Geneva) The b-quark mass including 0{\jm^) 

12:25 Lunch 
Afternoon Session, chair: Keh-Fei Liu 

14:30 Wolfgang Bietenholz Overlap hypercube quarks in the p- and the e-regime 

(Belin) 

15:05 Discussion: Lattice meets chiral perturbation theory * 

Hartmut Wittig (convenor) 
16:35 End of Session 
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Friday, October 6th, 2006 



Morning Session, chair: Christian Lang 



09:00 Keh-Fei Liu (Kentucky) 

09:35 Marina Dorati (Pavia) 
10:10 Andreas Schafer 
(Regensburg) 

10:20 

11:15 Volker Weinberg (Zeuthen) 
11:50 Pietro Faccioli (Trento) 
12:25 Christian Lang (Graz) 
12:30 



(j (600) as a tetraquark mesonium and 
the pattern of scalar mesons 

Moments of generalized parton distribution functions 
ChPT for generalized parton distributions (GPDs) 

Coffee 

The QCD vacuum as seen by overlap fermions 

The chiral regime of QCD in the instanton liquid model 

Farewell 

Lunch and End of Workshop 



no abstract was provided 

summary given at the end of these mini-proceedings 
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Results from Lattice QCD simulations with a twisted 

mass 

Luigi Scorzato 1 , on behalf of the ETM Collaboration 
X ECT*, strada delle tabarelle 286. 38050 - Villazzano (TN). Italy 

We report on first results of an ongoing effort to simulate lattice QCD with two 
degenerate flavours of quarks by means of the twisted mass formulation tuned to 
maximal twist. We obtain pseudo-scalar masses below 300 MeV on volumes with 
spatial size larger than 2 fm at values of the lattice spacing similar or smaller than 
0.1 fm. pp. We present first comparison with ChPT including Finite Size Effects. 
With the help of Wilson-ChPT we give evidence that 0(a 2 ) lattice artifacts are 
under control. Additionally, exploratory results for the case of Nf = 2 + 1 + 1 
flavours are discussed. 



References 



[1] K. Jansen and C. Urbach [ETM Collaboration] , |arXiv: hep-lat /0610015 
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Lattice QCD Simulations at Small Quark Masses 



Gerrit Schierholz 



Deutsches Elektronen-Synchrotron DESY 
D-22603 Hamburg, Germany 

John von Neumann-Institut fiir Computing NIC 
Deutsches Elektronen-Synchrotron DESY 
D-15738 Zeuthen, Germany 



Due to improvement in algorithms and computer performance unquenched 
simulations of Wilson-type fermions with lighter quark masses are now possi- 
ble P], [2]. This enables us to make contact with chiral perturbation theory 
(ChPT) and the real world. 

We report results for the pseudoscalar decay constants fp§ at pion masses 
down to 0(300) MeV, using Nf = 2 nonperturbatively 0(a) improved Wilson 
fermions. For calculational details see [3] . The renormalzation constants and im- 
provement coefficients of the axial vector current are computed nonperturbatively 
as well. 

We first investigate to see if we are entering a regime, where chiral logarithms 
are becoming visible, and look at the ratio of nondegenerate, partially quenched 
decay constants [4] 



(V: valence, S: sea quark). Our data displays chiral logarithms of about the 
expected size. Next we fit our data to the predictions of NLO ChPT. We ob- 
tain «4 m —0.58, «5 ~ —0.45. While is in reasonable agreement with other 
phenomenological estimates, is not. Using r = 0.5 fm [r = 0.467 fm] to set 
the scale, we find f K +/ 'f n + = 1.24 [fx+l fn+ = 1-22], to be compared with the 
experimental value 1.223, and /„+ = 76(4) (2) MeV = 81 (4) (2) MeV]. The 
first error is statistical, while the second error is due to the error in r /a. 



References 

[1] M. Gockeler et al, |hep-lat/0610066| 
[2] M. Gockeler et al, hep/lat/0610071. 

[3] M. Gockeler et al, Phys. Rev. D57 (1998) 5562 [hep/lat/9707021]. 

[4] S.R. Sharpe, Phys. Rev. D56 (1997) 7052, Erratum ibid. D62 (2000) 



For the QCDSF Collaboration 
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Hadrons for Chirally Improved Fermions 



T. Burch 1 , C. Gattringer 2 , L. Y. Glozman 2 , C. Hagen 1 , 
D. Hierl 1 , C. B. Lang 2 and Andreas Schafer 1 

(BGR [Bern-Graz-Regensburg] Collaboration) 
1 Inst. f. Theoretische Physik, Universitat Regensburg, D-93040 Regensburg, 
2 Inst. f. Physik, FB Theoretische Physik, Karl-Franzens-Universitat Graz, 

A-8010 Graz,Austria 

We present results for masses of excited mesons [I] and baryons [2] from a 
quenched calculation with Chirally Improved quarks at pion masses down to 
350 MeV. Our analysis of the correlators is based on the variational method. 
It was shown [3] that in this apporach ghost-channels may be identified safely. 
The key features are the use of a matrix of correlators from various source and 
sink operators and a basis which includes Jacobi smeared quark sources with 
different spatial widths, thereby improving overlap with states exhibiting radial 
excitations. In order to provide a large basis set for spanning the physical states, 
we also use interpolators with different Dirac structures. Our spectroscopy results 
for a wide range of ground state and excited hadrons are discussed. 



positive parity negative parity 




The figures show the results of straightforward extrapolations to the chiral region 
(mesons and positive as well as negative parity baryons). The horizontal bars 
represent the experimental numbers (where known). Filled symbols are used for 
those states where no corresponding state is listed in the particle data summary. 



References 

[1] T. Burch, C. Gattringer, L. Ya. Glozman, C. Hagen, C. B. Lang, and An- 
dreas Schafer, Phys. Rev. D 73 (2006) 094505, [arXiv: |hep-lat/0601026] . 

[2] T. Burch, C. Gattringer, L. Ya. Glozman, C. Hagen, D. Hierl, C. B. Lang, 
and A. Schafer, Phys. Rev. D 74 (2006) 014504, [arXiv: |hep-lat/0604019| . 

[3] Tommy Burch, Christof Gattringer, Leonid Ya. Glozman, Christian Hagen, 
and C. B. Lang, Phys. Rev. D73 (2006) 017502, [arXiv: |hep-lat/05 11054] . 
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Dynamical domain wall fermion simulations and the 
chiral perturbation theory 



Meifeng Lin 

Department of Physics, Columbia University, 
New York, NY 10027, USA 



The domain wall fermion formulation has the advantage of possessing exact 
flavor symmetry, and chiral symmetry is only mildly broken by a controllable 
amount. This allows us to compare the numerical results from the domain wall 
fermion simulations directly to the predictions of the continuum chiral pertur- 
bation theory [T], up to corrections from the residual chiral symmetry breaking, 
quantified as m res [2], and possible 0(a) lattice artifacts. From a field theoretic 
point of view, the residual mass m res serves as an additive mass renormalization 
to the input quark mass, and the chiral extrapolations are much simplified com- 
pared to the Wilson or Staggered fermion formulation. For domain wall fermions, 
the contamination from higher-dimension operators as a result of the 0(a) lattice 
artifacts is also exponentially suppressed by the span of the fifth dimension [3]. 
And the Symanzik effective action can be written as 



S eff = / d*x$(x)(i>y ll D' M - m q )^{x)\ + ae- aL °c dwf ^(x)a^F^(x) (1) 



where m q = nif + m rcs . Guided by this equation, I show that to the next-to- 
leading order, including the 0(a) chiral symmetry breaking term in the chiral 
expansions introduces only one new parameter [3]. 

I also present numerical results for the pseudoscalar masses (Mps) and de- 
cay constants {fps) from a series of 2+1 flavor dynamical domain wall fermion 
simulations carried out by the RBC and UKQCD collaborations [3]. When the 
Goldstone pion mass is above 400 MeV, the numerical results are not consistent 
with the NLO chiral perturbation theory in that simultaneous fits of Mp s and 
fps fail badly. However, for pion masses below 400 MeV, our preliminary data 
shows reasonable consistency with the predictions of NLO ChPT. 



References 

[1] J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465. 

[2] T. Blum et ai, Phys. Rev. D 69 (2004) 074502 [arXiv:hep-lat/0007038] . 

[3] M. Lin, PoS LAT2006 (2006) 185 |arXiv:hep-lat/061005"2 

[4] C. Allton et ai, PoS LAT2006 (2006) 096 
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Applications of Staggered Chiral Perturbation Theory 



C. Aubin 

Dept. of Physics, Columbia University, 
New York, NY USA 

One can account for the chiral symmetry violations that arise from staggered 
fermions on the lattice pp. Given current lattice simulations performed by the 
MILC Collaboration, the 0(a 2 ) taste-violating errors are of the same order as 
the light quark masses. Some results using staggered chiral perturbation theory 
(S^PT) when applied to light meson physics have been presented in Ref. [2]. 

This can be extended to include heavy quark effective theory with S%PT to 
understand the chiral behavior of lattice data for heavy- light quarks [3]. From 
this formalism, several results have been obtained in Ref. pE], most notably a 
prediction of the D decay constant before they were measured by CLEO-c plj. 

Applying these techniques to the kaon mixing parameter is also possible 
[6], however it becomes difficult to apply to lattice data due to the large number 
of free parameters. An option is to used a "mixed action" approach, where one 
uses domain- wall valence quarks and staggered sea quarks [7]. This technique is 
a slight modification of the continuum form for B% and is a promising method 
with which to extract this quantity. 



References 

[1] C. Bernard , Phys. Rev. D 65, 054031 (2002); C. Aubin and C. Bernard, 
Phys. Rev. D 68, 034014 (2003); Phys. Rev. D 68, 074011 (2003). 

[2] C. Aubin et al, Phys. Rev. D 70, 094505 (2004); C. Aubin et al. [Fermilab 
Lattice, MILC, and HPQCD Collaborations], Phys. Rev. D 70, 031504 
(2004); C. Aubin et al. [MILC Collaboration], P hys. Rev. D 70, 1145 01 
(2004); C. Bernard et al. [MILC Collaboration] , |arXiv:hep-lat/0609053~| 

[3] C. Aubin and C. Bernard, Phys. Rev. D 73, 014515 (2006), J. Laiho and 
R. S. Van de Water, Phys. Rev. D 73, 054501 (2006). 

[4] C. Aubin et al. [Fermilab Lattice and MILC Collaborations], Phys. Rev. 
Lett. 94, 011601 (2005); M. Okamoto et a/., Nucl. Phys. Proc. Suppl. 140, 
461 (2005); C. Aubin et al, Phys. Rev. Lett. 95, 122002 (2005). 

[5] M. Artuso et al. [CLEO Collaboration], Phys. Rev. Lett. 95, 251801 (2005). 

[6] Ruth S. Van de Water and S. R. Sharpe, Phys. Rev. D 73, 014003 (2006). 



[7] C. Aubin, J. Laiho and Ruth S. Van de Water, arXiv:hep-lat/0609009 



10 



Chiral Perturbation Theory and Lattice QCD with 

Wilson-type Quarks 

Sinya Aoki 

Graduate School of Pure and Applied Sciences, University of Tsukuba, 
Ten-oh-dai 1-1-1, Tsukuba, Ibaraki 305-8571, Japan 

I explain a necessity of introducing the lattice spacing effect into the chiral 
perturbation theory to fit lattice data obtained with Wilson-type quarks such as 
pion masses and decay constants. This formulation is called the Wilson chiral 
perturbation theory(WChPT) [T] . The WChPT has been applied to the analysis 
of the twisted-mass QCD [2] , [3] , [6] , [8] . The twisted- mass QCD becomes auto- 
matically 0(a) improved as long as the twisted-mass is set to the maximal value 
(the maximal twist). In actual numerical simulations, however, several different 
definitions of the maximal twist have been employed, and different definitions 
show quiet different chiral behaviours. I show that these differences among sev- 
eral definition of the maximal twist can be understood by the WChPT, and that 
pion mass and decay constant in quenched twisted-mass QCD can be fitted as 
a function of quark mass by the WChPT very well. The WChPT has been ex- 
tended to Nf = 2 + 1 flavor QCD for pseudo-scalar meson massesjl] and for 
vector meson masses [5]. These formula have been applied to Nf = 2 + 1 flavor 
QCD, in order to fit data. I show some examples of some chiral fits [7]. 



References 
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Aoki, Phys. Rev. D68 (2003) 054508 
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Aoki, O. Bar, Phys. Rev. D70 (2004) 116011 


[3] 
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Aoki, O. Bar, PoS LAT2005 (2005) 046 
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Aoki, O. Bar, T. Ishikawa, S. Takeda, Phys. Rev. D73 (2006) 014511 
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Aoki, O. Bar, S. Takeda, Phys. Rev. D73 (2006) 094501 
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Aoki, O. Bar, Phys. Rev. D74 (2006) 034511 
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Lattice QCD with mixed actions: 
Overlap fermions on a twisted mass sea 



Oliver Bar 1 , K. Jansen 2 , S. Schaefer 2 , L. Scorzato 3 , A. Shindler 2 

institute of Physics, Humboldt University, 
Newtonstrasse 15, 12489 Berlin, Germany 

2 NIC, DESY, 
Platanenallee 6, 15738 Zeuthen, Germany 

3 ECT* 

Strada delle Tabarelle 286, 38050 Villazzano (TN), Italy 

We present first results of a mixed action project [T]. We analyze gauge config- 
urations generated by the ETM collaboration (for recent summaries see [2], [3]). 
Two flavors of dynamical twisted mass fermions are taken into account with 
mps ~ 300 MeV. Neuberger's overlap Dirac operator [I] is used for the valence 
sector. For the quark mass matching we tune the charged valence pion mass such 
that it equals the charged sea pion mass. 

In order to check the partial quenching effects in our mixed action theory we 
compute the scalar correlator. As expected, the correlator is negative [5] . Mixed 
action theories can also be studied with ChPT methods [S] , [Z] , [H] • We fit our 
data to the LO ChPT result for the scalar correlator of Ref. [9] and find good 
qualitative agreement. 
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Thoughts on chiral extrapolations for excited hadrons 

Ulf-G. Meifiner 1 * 2 

1 HISKP (Th), Universitat Bonn, D-53115 Bonn, Germany 
2 IKP (Th), Forschungszentrum Jiilich, D-52425 Jiilich, Germany 

An understanding of QCD without ist excitation spectrum is incomplete. In 
this talk I review work on chiral extrapolations for excited hadron states - I con- 
centrate on the p(770) and the Roper A/"* (1440). Both these fields are considered 
as explicit degrees of freedom in the pertinent chiral effective Lagrangian. The 
new (large) scale related to their mass introduces complications in the power 
counting, how to treat this in case of the real part of the rho-meson self-energy is 
discussed in [I] and the even more complicated nucleon-Roper-pion system is con- 
sidered in [2] by extending the standard infrared regularization scale. The quark 
mass expansion of the p looks very similar to the one of the nucleon, analyzing 
e.g. the old CP-PACS data, we find 650 MeV < M° < 800 MeV for the chiral 
limit mass. For the Roper, the quark mass expansion comes out very close to the 
one of the nucleon. With coupling constants of natural size, no sharp decrease 
of mR oper (M 7r ) for small pion masses is observed. Of course, these considerations 
need to be sharpened. For the p, one has to analyze the more recent data and in 
case of the Roper, an extension to include the Air and N(tttc)s channels is called 
for. The A(1232) resonance is discussed in Bernard's talk [3] and for dynamically 
resonances I refer to Oiler's talk [3]. 
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Chiral extrapolation of baryon properties 

V. Bernard 

Laboratoire de Physique Theorique, 67085 Strasbourg, France 

Enormous progresses have been made in lattice calculation in the last years. 
However these calculations still involve rather large pion masses. It is thus nec- 
essary to perform chiral extrapolations in order to make contact to the physical 
world. This can be done in a model independent manner using chiral perturba- 
tion theory. For discussion on the different regularization used in baryon CHPT 
see [1] . In this talk I reported on chiral extrapolation of some baryon properties, 
namely the mass [2] and axial vector coupling [3] of the nucleon and the A mass 
[I] essentially in the continuum limit. In the nucleon case two loop calculations 
are by now available (only the 0(p 5 ) terms have been determined in the case of 
the mass). In these calculations only 7 combinations of LEC's from the ttN La- 
grangian up to dimension 4 appear, 5 of these being rather well determined from 
ttN scattering and the Goldberger Treiman discrepancy. Taking into account the 
error bars on these LECs and studying the convergence of the series one finds 
that the chiral extrapolation of these two quantities can be trusted up to ~ 350 
MeV. Taking explicitely into account the A degree of freedom does not improve 
on the result [5]. Discussion on finite volume effects can be found for exemple in 
[6]. In the case of the A mass fits to the lattice data in the small scale expansion 
seems to indicate that the symmetry breaker LEC % is far from its SU(6) value 
leading to a irA sigma term smaller than the ttN one. Further study is needed to 
confirm this result. The conclusion of these studies is that it is not usefull to fit 
lattice data up to rather large pion masses but that it is much more important 
to perform simultaneous systematic chiral extrapolation with a consistent set of 
LEC's and to carefully evaluate the theoretical uncertainties. 
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Some topics in manifestly Lorentz- invariant ChPT 

D. Djukanovic 1 , J. Gegelia 1 , M. R. Schindler 1 and S. Scherer 1 

1 Institut fur Kernphysik, Johannes Gutenberg-Universitat 
55099 Mainz, Germany 

We present three topics in manifestly Lorentz-invariant baryon chiral pertur- 
bation theory (BChPT). 

The higher-derivative formulation of ChPT pQ presents an alternative regu- 
larization scheme that can be used in the vacuum, one- and few nucleon sector 
of ChPT. By performing field transformations in the canonical Lagrangian, ad- 
ditional higher-derivative terms are introduced in the lowest-order Lagrangian. 
These terms improve the ultraviolet behavior of propagators while preserving all 
symmetries. The new parameters can be interpreted as smooth cutoffs, whose 
values can be chosen freely. 

To investigate the applicability of chiral expansions [2] we compare the sum 
of an infinite number of terms contributing to the nucleon self energy starting at 
the two- loop level with the leading non- analytic term ~ M 3 . For pion masses 
M > 550 MeV the contributions from the sum of higher-order terms is larger 
than the leading-order contribution. This shows that for M > 550 MeV the 
power counting is no longer valid. 

We also present the results of the first complete two-loop calculation of the 
nucleon mass [3] using the reformulated infrared renormalization [4] and compare 
our result to the heavy baryon result at order 0(q 5 ) [5]. While at order 0(q 5 ) only 
low energy constants appear that have previously been determined, the terms at 
order 0(q 6 ) contain combinations of numerically unknown constants. 
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Chiral Extrapolations in Few— Nucleon Systems 

Evgeny Epelbaum 1,2 
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The dependence of the low-energy dynamics of QCD on the variation of the 
up and down quark masses can be naturally studied within the framework of 
chiral effective field theory. While the situation in the Goldstone boson and 
single-nucleon sectors is rather straightforward due to the validity of perturbation 
theory, it is highly non-trivial in the few-nucleon sector where non-perturbative 
methods are required to understand e.g. the properties of the shallow bound 
states. The quark mass dependence of various few-nucleon observables is not 
only of academic interest, but also relevant for interpolating the results from 
lattice gauge theory [I], [2] and imposing bounds on the time-dependence of 
fundamental couplings. 

We have studied the quark mass dependence of the nucleon-nucleon (NN) 
force and various two-nucleon observables at next-to-leading order in the chiral 
expansion [3], see also [I] for a similar work. We found that both NN S-wave 
scattering lengths become smaller in magnitude (i.e. more natural) in the chiral 
limit. The deuteron is less bound when the quark mass is increased and becomes 
unbound for M v larger than ~ 200 MeV. The theoretical uncertainty of our results 
is rather large and mainly caused by the uncertainty in the determination of the 
LEC g? 16 and by the unknown LECs associated with the leading M^-dependent 
NN contact interactions. We have also studied the infrared renormalization group 
limit cycle in the three-nucleon system [5], whose existence for certain values of 
m u and md was conjectured in [6]. 
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Nuclear Lattice Simulations 

B. Borasoy 

Helmholtz-Institut fur Strahlen- und Kernphysik (Theorie) 
Universitat Bonn, NuBallee 14-16, D-53115 Bonn, Germany 

Nuclear lattice simulations are a novel approach to nuclear physics based 
on numerical simulations of chiral effective field theory on the lattice. Within 
this scheme the chiral effective Lagrangian is formulated on the lattice and the 
corresponding path integral is evaluated numerically by Monte Carlo sampling. 

In two exploratory studies we have employed this framework to calculate prop- 
erties of the lightest nuclei: the deuteron, triton and Helium-4. In [1] the triton 
was investigated in the SU(4) limit of the pionless theory including the three- 
nucleon force. By comparing with continuum results we demonstrated that the 
nuclear lattice formalism can be used to study few-body nucleon physics. A first 
measurement of the three-nucleon force on the lattice in the Wigner symmetry 
limit has been provided. Moreover, we have shown that many-body simulations 
of cold dilute nuclear matter in the Wigner symmetry limit should be possible 
without a sign problem. This work is thus also of relevance for future many-body 
simulations with arbitrary numbers of nucleons including three-body effects. 

The pionful theory is considered in [2] where we calulate the binding energies 
and sizes of triton and 4 He. A Monte Carlo algorithm with pions and auxiliary 
fields is constructed which reproduces both the lowest order S'-wave contact in- 
teractions and instantaneous one-pion exchange between the nucleons. Due to 
the approximate Wigner symmetry of the actual interactions we find only mi- 
nor sign and phase oscillations for light nuclei. The importance of higher order 
interactions is also investigated. 

Although our investigations are still at the exploratory level, lattice simula- 
tions of chiral effective field theory appear to be a promising tool to investigate 
few- and many-body nuclear physics with a clear theoretical connection to QCD 
and a systematic chiral expansion. 
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Chiral Non-Perturbative Study of Pseudoscalar 

Self-Energies 

Luis Roca and Jose Antonio Oiler 

Departamento de Ffsica, Universidad de Murcia, 
E-30071 Murcia, Espana 

In this talk we report about our recent work [I]. We perform a non-perturbative 
chiral study of the masses of the lightest pseudoscalar mesons. The pseudoscalar 
self-energies are calculated by the evaluation of the scalar self-energy loops with 
full S-wave meson-meson amplitudes taken from Unitary Chiral Perturbation 
Theory (UCHPT) [2]. These amplitudes, among other features, contain the 
lightest nonet of scalar resonances a, /o(980), a (980) and k and the heavier 
ones up to 1.5 GeV. The self-energy loops are regularized by a proper subtrac- 
tion of the infinities within a dispersion relation formulation of the scattering 
amplitudes. Values for the bare masses of pions and kaons are obtained as well 
as an estimate of the mass of the r/g. We then match to the self-energies from 
standard Chiral Perturbation Theory (CHPT) to 0(p 4 ) and resum higher orders 
from our calculated scalar self-energies. The dependence of the self-energies on 
the quark masses allows a determination of the ratio of the strange quark mass 
upon the mean of the lightest quark masses, m s /m, in terms of the 0(p A ) CHPT 
low energy constant combinations 2Lg — L\ and 2Lg — L\. In this way, we give 
a range for the values of these low energy counterterms and for 3L7 + Lg, once 
the 77 meson mass is invoked. The low energy constants are further constraint by 
performing a fit to the recent MILC lattice data on the pseudoscalar masses. An 
excellent reproduction of the MILC data is obtained, at the level of 1% of relative 
error in the pseudoscalar masses, and m s /rh = 25.6 ± 2.5 results. This value is 
consistent with 24.4 ±1.5 from CHPT and phenomenology and more marginally 
with the value 27.4 ± 0.5 obtained from pure perturbative chiral extrapolations 
of the MILC lattice data to physical values of the lightest quark masses. We also 
show that our estimation of the 0(p 6 ) and higher chiral orders, indicates that 
the SU(3) CHPT series on pseudoscalar self-energies seems to stabilize and be- 
have much better once the 0(p 6 ) CHPT contributions to the pseudoscalar masses 
are included, despite that the 0(p 4 ) contributions are typically smaller than the 
0(p 6 ) ones. 
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Scalar two-loop diagrams on the lattice 

Hermann Krebs 

HISKP, Universitat Bonn and IKP, Forschungszentrum Jiilich 

Numerical simulations on the lattice are always performed in finite volume 
with a finite lattice spacing. In order to understand both the infinite volume and 
the continuum limit non-trivial extrapolations are needed. In some cases lattice 
perturbation theory provides a powerful tool in getting an analytical control over 
nonperturbative simulations. 

Lattice perturbation theory itself is a challenging field. Since Lorentz-invariance 
is lost on the lattice many tools developed for continuum perturbation theory 
such as Feynman parameters are not applicable on the lattice. Especially for 
multi-loop calculations where conventional methods often lead to loss in preci- 
sion new techniques are highly desirable. A novel powerful technique in this 
respect was introduced by Luscher and Weisz [T]. Instead of integration over 
the first Brillouin-zone in momentum space they suggested to calculate infinite 
sums in coordinate space. In this case extensive knowledge of the free lattice 
propagators is required. Recursion relations for the massless propagator were 
derived which allowed to express the Green function at any lattice site as a linear 
combination of two constants which can be determined to very high precision. 

We extended this technique to the massive case and could give the small 
mass expansion of the massive scalar propagator [2]. In the BPHZ subtraction 
scheme any scalar two loop diagram can be expressed as a sum of a regular and 
a singular part which converge and diverge in the continuum limit, respectively. 
The singular part is always of the following sunset type: 

21 

J-n (2VT) 4 (27r) 4 p + m 2]a[£2 + ^[(f^q) + m 2] 7 ' '' V 2 

We developed a method for deriving the small mass expansion of any singular 
part [3]. The coefficients of this expansion can be calculated numerically to very 
high precision such that even discretization effects can also be studied within our 
method. The proposed method for the small mass expansion is universal and can 
be applied once the particle propagator is known, e.g., it can easily be extended 
to Wilson fermions. 

References 

[1] M. Luscher, P. Weisz, Nucl. Phys. B 445 (1995) 429. 
[2] B. Borasoy, H. Krebs, Phys. Rev. D 72 (2005) 056003. 
[3] B. Borasoy, H. Krebs, Nucl. Phys. B 748 (2006) 1. 



19 



Nucleon and N to A form factors in Lattice QCD 

C. Alexandrou 1 , Th. Leontiou 1 , J. W. Negele 2 and A. Tsapalis 3 
x Dep. of Physics, University of Cyprus, CY-1678 Nicosia, Cyprus 
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The isovector nucleon form factors are evaluated in lattice QCD in the quenched 
theory and using two degenerate flavors of dynamical Wilson fermions for a range 
of pion masses between 690 and 380 MeV, using the nucleon mass to set the lattice 
spacing a. We find that the isovector electric form factor, Ge, at the physical limit 
deviates more from experiment than the magnetic form factor, Gm [I]- Lattice 
results show a weaker g 2 -dependence for the isovector ratio (j,Ge/Gm as compared 
to the results obtained in recent polarization experiments [2]. The electromag- 
netic and axial N to A transition form factors are evaluated for quenched and 
two dynamical flavor of Wilson fermions as well as in a hybrid scheme where we 
use MILC configurations and domain wall fermions over a range of pion masses 
between 690 and 360 MeV. We find that lattice results for the dominant magnetic 
dipole form factor in 7*iV — > A, when linearly extrapolated in m 2 to the physical 
point, yield larger values than experiment [3]. Results for the ratios of electric 
and coulomb quadrupole amplitudes to the magnetic dipole, EMR and CMR, are 
non-zero and negative as in experiment albeit with statistical errors that are still 
too large in the unquenched case [I] to allow assessment of pion cloud contribu- 
tions. Our results for the four N to A axial form factors show that C% ~ 0, 
is small and and C$ are the dominant form factors. We evaluate the ratio 
which is the analog of gA/gv, as a function of Q 2 . This ratio provides 
the leading contribution to the the parity violating asymmetry [5] to be measured 
by the GO experiment at JLab [6]. The off-diagonal Goldberger-Treiman relation 
is examined by evaluating, at the same lattice parameters, g n NA and f^. 
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Hadron Structure using DWF Quarks on an Asqtad Sea 

David Richards, for Lattice Hadron Physics Collaboration 
Jefferson Laboratory, 12000 Jefferson Avenue, Newport News, VA 23606, USA. 

Moments of unpolarized, helicity, and transversity distributions, electromag- 
netic form factors, and generalized form factors of the nucleon are presented from 
a preliminary analysis of lattice results using pion masses down to 359 MeV[T]. 
We employ a hybrid approach, in which improved, staggered quarks are used 
for the generation of the gauge configurations, whilst domain-wall fermions, with 
their desirable chiral properties, are used for the valence quarks. 

The nucleon axial-vector charge, a benchmark quantity of QCD, is particularly 
robust under chiral extrapolation; the consistency of the hybrid calculation, both 
with other lattice calculations, and with experiment at the physical pion mass, 
is encouraging [2]. Lattice moments of structure functions and GPDs likewise 
require extrapolation to the physical quark masses; a long-standing puzzle has 
been the flat behaviour of the flavour- non-singlet momentum fraction, (x), of 
the nucleon, at a value considerably higher than the experimental value. An 
approach in which we apply xPT, with low-energy constants g& and given by 
their lattice values at each quark mass, allows a two-parameter extrapolation in 
m 1 ^ I to yield a value for (x), and other benchmark quantities, at the physical 
quark masses that are consistent with experiment. This development encourages 
to now exploit the predictive power of these calculations. 

The low-Q 2 behaviour of the nucleon form factors describes the distribution 
of charge and magnetism within a nucleon. The slope of the Fi form factor is 
related to the rms charge radius; the chiral extrapolation of the isovector charge 
radius likewise yields values consistent with experiment [3] . Generalized Parton 
Distributions provide new insight to hadron structure. For example, the total 
angular momentum carried by the quarks is related to a combination of moments 
J q = 1(^20^+^20^) EO- Combined with measurements of quark spins, we find the 
total orbital angular momentum carried by quarks is small, though that carried 
by individual flavours is substantial. 
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Distribution amplitudes from the lattice 
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Hadronic wave functions are of crucial importance when describing exclu- 
sive and semi-exclusive reactions [I]. For detailed references and applications 
consult [2J. Distribution amplitudes (DAs) are related to the hadron's Bethe- 
Salpeter wave function, <j)(x,k±), by an integral over transverse momenta. For 
the leading twist meson-DAs we have 

0(x,/i 2 ) = Z 2 (/j 2 ) / d 2 k(f)(x,k ± ) , (1) 

J |fcx|<M 

where x is the quark longitudinal momentum fraction, Z 2 the renormalization fac- 
tor (in the light-cone gauge) for the quark-field operators in the wave-function, 
and n denotes the renormalization scale. In this presentation we quote all num- 
bers with a scale /x 2 = 4 GeV 2 in the MS-scheme. 

It is convenient to rescale £ = 2x — 1. It is common to expand DAs into 
their Gegenbauer moments and quote the expansion coefficients, at a given 
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renormalization scale parameterization of DAs, 

<Kt, = |(1 - e 2 ) ^l + g ^(/^(o) ■ (2) 

The zeroth moment is normalized to unity, J_j c?£0(£, ji 2 ) = 1, at any energy scale 
/i 2 . Taking the u- and <i-quarks to be degenerate, G-parity implies that the pion 
DA is an even function of £ and hence all odd moments vanish, i.e., a^n+i = 0- 

Recently, we have computed the first moments of meson distribution ampli- 
tudes [2] on the lattice. Independently, a calculation of the first moment of the 
kaon distribution amplitude has appeared which uses a different discretization 
scheme and different working points [3]. We find that their results are compati- 
ble with ours. 

This presentation discusses the calculation of the first non- vanishing moment 
of the pion DA, aj, and the first two moments of the kaon DA, af and ■ We 
compare our results to previous estimates from sum rules and experiment and 
discuss the phenomenological implications of our lattice computation. 

We obtain <22(4GeV 2 ) = 0.201(114) which disfavors the commonly adopted 
"asymptotic" model with a\ n = 0. This result is also compatible with model 
estimates, cf. e.g. [I]. 

We also find af-(4GeV 2 ) = 0.175(18)(47) which provides an important test 
of different and competing kaon models employed in the literature, see [5]. 

Finally, we find a^(A GeV 2 ) = 0.0453(9) (29) which not only confirms the sign, 
but also the order of magnitude for the asymmetry of the kaon wave function. 
This key finding is of major importance both for the understanding of B physics 
and of weak decays [5]. 
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Pion structure from lattice QCD 
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We calculate moments of generalised parton distributions (GPDs) of the pion. 
This is done using a simulation of non-perturbatively improved dynamical Wil- 
son fermions with two flavours with Wilson glue. The configurations have been 
generated within the QCDSF/UKQCD/DIK collaborations and cover a range of 
pion masses down to ~ 340 MeV. We show partly preliminary results for the 
lowest (two) moment(s) of the tensor (vector) GPD H^ n {H%), defined as 

2 pa* /"HA A A 

— H*= / —e iXn - Px (n(p>)\q(--n)rUqhn)\n(p)) + higher twist, (1) 
n ■ r J Z7T 2 2 

^7p H h= y ^e iA "^(7r(pO|g(-^)i^Wg(-n)|7r(p))+higher twist. (2) 

The moments of these GPDs are parametrised by generalised form factors. We 
focus on the lowest moment of the vector GPD Eq. (TjQ) which is the pion form 
factor Details of the calculation and the lattice setup and further references 
can be found in P],[2]. The pion form factor can be very well described by a 
monopole ansatz. The squared monopole masses obtained from this ansatz are 
extrapolated to the physical point linearly against square of the pion mass. From 
that we obtain a charge radius of (r 2 ) = 0.440(19) fm 2 in good agreement with 
experiment. We also perform a study of the finite size and scaling effects [2]. 
The second moment of the vector GPD corresponds to the momentum fraction 
of the quarks inside the pion. A preliminary result indicates that quarks and anti- 
quarks inside the pion carry ~ 50% of the momentum fraction. The calculation 
of the pion tensor GPD Eq. (J2J) has not been attempted before. This quantity is 
experimentally unknown and thus provides genuine new insights to the internal 
spin structure of the pion [3] . Our preliminary data so far indicate a rich structure 
with a no n- vanishing tensor GPD also in the forward limit. 
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Nucleon Structure Functions and Form Factors 

Dirk Pleiter 1 
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In this talk we have reported on recent results for the moments of the unpo- 
larised nucleon structure functions pQ and the electromagnetic form factors [2J. 
The calculations have been done using dynamical gauge field configurations with 
Nf — 2 flavours of non-perturbatively 0(a)-improved Wilson fermions, which 
have been generated by the QCDSF, DIK and UKQCD collaborations. On the 
different gauge field ensembles the pseudo-scalar meson mass is in the range of 
340 to 1170 MeV. The lattice spacings vary between 0.07 and 0.11 fm, i.e. our 
lattice are relatively fine. The spatial extension of our lattice is 1.5 fm for the 
heavy sea quark masses and up to 2.6 fm for the very light quark mass. The 
renormalisation constants have been determined non-perturbatively. 

By comparing our results for different lattice spacings and different volumes 
we find that discretisation and finite volume effects are relatively small. The 
systematic error which seems most difficult to control stems from the extrapo- 
lation of the lattice results to the physical mass of the up/down quarks. When 
using a naive ansatz linear in the quark mass, which describes our data actually 
very well, we find a significant discrepancy with the experimental values. Using 
the quenched approximation we previously obtained surprisingly similar results 

For the moments of the structure functions as well as for the form factor 
radii and anomalous magnetic moments chiral effective theories have been used 
to investigate the quark mass dependence. These calculations suggest a rather 
strong quark mass dependence at very small quark masses. The relevant range 
of quark masses is only starting to become accessible for lattice simulations. 
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The kaon bag-parameter, B K = {K°\s^^l-^)ds^^l-^)d\K°) / (y^m^), 
defined as the ratio of the AS* = 2 matrix element to its value obtained in the 
vacuum saturation approximation is a measure of indirect CP- violation and pro- 
vides an important non-perturbative input to the Cabbibo-Kobayashi-Maskawa 
matrix fit, cf. e.g. pQ. 

We present a preliminary study of the kaon bag-parameter Bk measured on 
L 3 x T x L s = 16 3 x 32 x 16 lattices using iVy = 2 + 1 dynamical flavors of 
domain wall fermions, see also [2], [3], [I]. The mass of the heavier (single) 
flavor was tuned to match the physical strange quark mass, while three different 
values for the masses of the two degenerate light quarks have been used in the 
simulation. These correspond to pion masses in the range of 390 to 630 MeV. 
Whereas the applicability of chiral perturbation theory to extrapolate the pion 
masses and decay constants to the physical limit is at least questionable in this 
mass range (see [5], [6] for a discussion), we still may use the partially quenched 
chiral perturbation theory formulae given in [7] as a smooth fitting function for 
an extrapolation of our results concerning Bk- Ignoring small non-degenerate 
effects even allows for an interpolation over just a small mass range. For details 
we refer to [3] and our upcoming paper [I]. 
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Low- energy constants for AS = 1 transitions 
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In order to investigate the origins of the AI = 1/2 rule in K — > tttt de- 
cays we have devised a strategy to determine the LECs of operators describing 
AS = 1 transitions pQ. We keep an active charm quark, express K — > tttt am- 
plitudes in terms of LECs and study their dependence on m c . By employing 
Neuberger fermions, the severe problem of mixing with lower dimensional oper- 
ators is completely avoided. The LECs are extracted by matching lattice QCD 
data for suitable ratios of correlation functions to the expressions of ChPT in 
finite volume [2J. The first step in our implementation of the strategy is focused 
on the case of a light and degenerate charm quark, i.e. the GIM limit (m c = m u ). 
In order to tame the extreme statistical fluctuations encountered in simulations 
near the chiral limit, noise reduction techniques (low- mode averaging) must be 
employed [3j. After renormalising ratios of correlators non-perturbatively [I], we 
perform joint fits to the corresponding ChPT expressions in both the e and p- 
regimes. Our results in the GIM limit [5] indicate a previously unseen, significant 
enhancement of the AI = 1/2 amplitude, which however, is still a factor 4 smaller 
than the experimental number. First attempts to analyse the direct influence of 
the charm quark in ChPT [5] will be supplemented by future numerical studies. 
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The aim of the investigation, carried out in the present work, is to study the 
feasibility of the extraction of the A-resonance parameters from the lattice data 
at the physical value of the pion mass, i.e. in the case of unstable A. To this 
end, we calculate the self-energy of the A-resonance in a finite volume up to and 
including 0(p 3 ) in Chiral Perturbation Theory, using infrared regularization. The 
poles of the A-propagator, obtained in these calculations, determine the energy 
spectrum of the correlator of two A-fields in a finite Euclidian box and can be 
obtained in lattice QCD. The results of our calculations enable one to investigate 
the dependence of the energy levels both on the size L of the box and on the pion 
mass M n . 

As previously argued in the literature (see, e.g. P,[2],[3]), the presence of 
a narrow resonance in the particle spectrum manifests itself through a peculiar 
irregular behavior of the energy levels with respect to L (avoided level crossing 
near the resonance energy). We have however verified by explicit calculations 
at 0(p 3 ) that, in the case of interest, the dependence of energy levels is rather 
smooth - the avoided level crossing is nearly washed out due to a large A-width. 
Despite this fact, on the basis of the numerical study, we argue that the extraction 
of the parameters of A (namely, of the mass and the irNA coupling constant) 
from a measured L-dependence of the lowest energy levels is still a feasible task. 
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The QCDSF collaboration (together with the DESY-ITEP-Kanazawa (DIK) 
collaboration) is generating gauge field configurations with two degenerate flavours 
of nonperturbatively 0(a)-improved Wilson fermions (clover fermions) and Wil- 
son glue, including in the analysis also older ensembles from UKQCD. In this 
talk we have discussed our approach to the treatment of finite size effects. It is 
based on appropriate versions of chiral effective field theory and hence intimately 
connected with the quark mass dependence of the physical quantities considered. 
In particular we deal with data from our two finite size runs (three volumes each, 
with the other simulation parameters kept fixed). Ongoing simulations at smaller 
quark masses are expected to shed new light on the chiral behaviour (see, e.g., 
Refs. P,[2]). 

Among the quantities studied are the nucleon mass [3] and the axial coupling 
constant of the nucleon [4]. For the pion (or pseudoscalar) mass and the pion 
decay constant we make use of the recently proposed "resummed" Luscher type 
formulae [5]. 
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In this talk we explore the possibilities of chiral extrapolations using the covariant 
formulation of chiral effective field theories (EFTs) with baryons. In particular, a 
study of the convergence of the chiral expansion and of the statistical uncertain- 
ties arising from presently unknown parameters in the extrapolation functions 
is performed for the examples of the nucleon mass and the anomalous magnetic 
moment of the nucleon. In the calculation of these quantities at next to leading 
one loop level, i.e. order p 4 , we rely on a modified scheme of infrared regulariza- 
tion (IR), which will be discussed in-depth in a forthcoming publication. 
Calculating an observable in EFT, one has a freedom of choice which part of 
the analytic pieces of the result is considered to arise from loop dynamics and 
which one corresponds to local operators. Taking advance of this freedom, the 
IR-scheme was introduced [1] to overcome the long-known problem of power- 
counting violations in the MS-scheme of covariant BChPT. Compared to stan- 
dard MS-results, IR absorbs an infinite sum of terms analytic in the quark mass 
into low energy constants (LECs). Problem^ can arise due to the fact, that in 
EFT at a certain order only a finite number of LECs is included. Therefore, we 
propose to only consider those analytic terms as part of short distance physics, 
which can be absorbed into LECs present at the order of the calculation. 
The quark mass dependence both for the nucleon mass and for the anomalous 
magnetic moment of the nucleon calculated with this method has the following 
properties: 1. A successful chiral extrapolation of currently available lattice data 
to the physical point can be performed with reasonable values for the appearing 
parameters, i.e. values consistent with known information from scattering theory. 
2. Statistical uncertainties are small. 3. We do see a clear convergence of the 
results when going from p 3 to p 4 . 4. Finally, we estimate possible effects arising 
at next to leading order by using realistic quark mass dependencies for the lowest 
order coupling^ and find only little impact on the best-fit curves. 
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Chiral perturbation theory determines the low-energy structure of Green func- 
tions in QCD through a systematic expansion in powers of quark masses and of 
the external momenta [3],[2]- The method has been successfully applied for an 
expansion in powers of m u and at fixed m s ,m c , . . . [chiral SU(2) xSU(2)] as 
well as for an expansion in powers of m u , and m s at fixed m c ,mb,m t [chiral 
SU(3)xSU(3)]. It makes use of an effective action who's pertinent low-energy 
constants (LECs) encode the heavy degrees of freedom that are not explicitly 
contained in it. The LECs in SU(2) xSU(2) depend on the strange quark mass, 
while the ones in SU(3)xSU(3) are independent thereof. 

If one limits the external momenta to values small compared to m s and treats 
m u and as small in comparison to m s , the degrees of freedom of the K- 
and rj- mesons freeze and one can work out explicitly the strange quark mass 
dependence of the LECs in SU(2)xSU(2) from the effective action pertaining to 
SU(3)xSU(3). For example, for the pion decay constant F at m u = = 0, 
m s ^ 0, one has at one loop [2] 

+ • 

Here, Fo,Bq and L\ are LECs from SU(3) xSU(3), and the renormalization scale 
is denoted by \i. Relations of this type generate constraints on the possible 
values of the low-energy constants: information on LECs in SU(2) xSU(2) may 
be translated into information on LECs in SU(3) xSU(3) and vice versa [2J. 
In my talk, I discussed recent efforts [3] to establish the corresponding relations 
between the LECs to two-loop accuracy. 
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A Chiral Two-Matrix Theory and the Chiral Lagrangian 
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In the e-regime of QCD an efficient method for extracting the chiral conden- 
sate E is based on the distributions of Dirac operator eigenvalues. Previously, 
the analogous technique for the pion decay constant F n was based on introducing 
a vector source ("imaginary isospin chemical potential"). Compact analytical 
expressions have been found for the quenched case pQ and that of two light dy- 
namical quarks [2]. It is computationally difficult to go beyond these two cases. 

One would like to find a Random Matrix Theory [3] that describes the inclu- 
sion of imaginary isospin chemical potential. Because it seeks the microscopic 
spectra of Dirac operators with different chemical potential for u and d type 
quarks, it turns out to be given by a Random Two-Matrix Theory [4]. We solve 
this theory for finite size N of the matrices by means of biorthogonal polynomials. 
By deriving the relevant kernels we find explicit formulas for all (n, k) mixed or 
unmixed spectral correlation functions of the two involved Dirac operators. In the 
microscopic scaling limit we compute the corresponding analytical expressions as 
well. These correspond to the leading expressions of the e-expansion based on the 
chiral Lagrangian. In the two cases where comparison to the chiral Lagrangian 
framework are possible pQ, [2] we find exact agreement. A useful by-product 
is the corresponding analytical expressions for spectral correlation functions of 
Dirac operators with chemical potential, evaluated in ensembles with dynamical 
quarks of zero chemical potential. These expressions allow for a determination of 
the pion decay constant by means of ordinary lattice configurations. 
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The CKM element is important for the phenomenology of flavor physics in 
determining the apex of the unitarity triangle in the complex plane. It is possible 
to determine \V c b\ from both inclusive and exclusive semileptonic B decays, and 
they are both limited by theoretical uncertainties. The inclusive method makes 
use of the heavy quark expansion PQ, [2], but is limited by the breakdown of local 
quark-hadron duality, the errors of which are difficult to estimate. The exclusive 
method requires reducing the uncertainty of the form factor Tb^d*^ which has 
been calculated with lattice QCD in the quenched approximation using the double 
ratio method [3]. Given the phenomenological importance of this quantity we 
have revisited this calculation of Tb^d* using the 2+1 flavor MILC lattices with 
improved light staggered quarks @]. The lattice calculation was done on the 
MILC coarse lattices (a ~ 0.125 fm) where the light quarks were computed with 
the "asqtad" action. The heavy quarks were computed using the clover action 
with the Fermilab interpretation in terms of HQET [5] . We have looked at three 
light masses at the full QCD points, m va i ence = m sea . The lightest such mass was 
around m s /7. The chiral perturbation theory (ChPT) for heavy light mesons with 
a Wilson-like heavy quark and a staggered light quark was worked out by Aubin 
and Bernard in [6], and the ChPT relevant for B — > D* was done in [7]. The plan 
for the future is to add statistics by running on multiple time sources, and to run 
on additional lattice spacings generated by MILC and the Fermilab lattice group 
in order to determine the lattice spacing dependence. We will incorporate the 
matching coefficients into the analysis, and of course, we will do a careful analysis 
of all the systematic errors in order to report a final number for the B — ► D*tv 
form factor. 
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The b-quark is too heavy to be treated dynamically on the lattice. Indeed, 
in units of nowadays affordable resolutions am^ » 1. A theoretically attractive 
option is to use effective theories like Heavy Quark Effective Theory (HQET). 
Introduced in p], it provides the correct asymptotic description of QCD correla- 
tion functions in the limit — > oo. Sub leading effects are described by higher 
dimensional operators whose coupling constants are formally 0(l/rab) to the 
appropriate power. The degrees of freedom in the effective theory are strongly 
coupled and therefore a non-perturbative approach is needed. In addition HQET 
has to be matched to the full theory, a step, which must as well be performed 
non-perturbatively beyond leading order in l/m^. 

A framework for non-perturbative HQET on the lattice has been introduced 
in [2] . We present here a specific application of that approach to the computation 
of the b-quark mass in the quenched approximation. The details can be found 
in [3]. In a first step we match HQET and QCD in a small volume (L = 0.4 fm). 
To make contact with phenomenology the HQET expressions of the relevant 
quantities are then evolved to large volume. The essential tools in this step are 
the step scaling functions, which we introduce to describe the effects of a change 
in the linear size L of the system by a factor two. Eventually the b-quark mass 
is fixed using the (spin averaged) B s meson mass. We include the subleading, 
1/mb, terms in each step. Our final result is m^m^) = 4.347(48) GeV in the 
MS scheme. We have implemented twelve different matching conditions for this 
computation. Although the contributions at a given order may differ in the 
twelve different cases, once the LO and the NLO terms are summed together all 
the estimates agree. That implicitly suggests that the I /ml corrections to the 
quark mass are extremely small as expected. We plan to extend the method to 
the computation of the decay constant Fq s , which we have already computed in 
the static approximation [4] and to real (unquenched) QCD. 
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The overlap hypercube fermion is a variant of a chirally symmetric lattice 
fermion, which has a much higher degree of locality than the standard overlap 
fermion. We applied this formulation in quenched QCD simulations with light 
quarks [I], which were performed with HLRN machines on a 12 3 x 24 lattice 
at (5 — 5.85 (corresponding to a lattice spacing of a ~ 0.123 fm according to 
the Sommer scale). The topological susceptibility that we obtained supports the 
Witten-Veneziano scenario to explain the mass of rj'. 

In the p-regime we evaluated the masses of light pseudoscalar and vector 
mesons. Although we reach m n < 300 MeV it is difficult to extrapolate the 
corresponding value for F n — measured from a matrix element — to the physical 
point. However, we found a stable axial current renormalisation constant 
close to 1, in contrast to the results with the standard overlap operator. 

The use of tiny bare quark masses in the same volume took us to the e- 
regime, where we evaluated the leading Low Energy Constants £ and F n with 
various methods [1] . The densities of low lying Dirac eigenvalues were fitted to 
predictions by chiral Random Matrix Theory, which yields X ~ (300 MeV) 3 . 
Next we compared the axial current correlators to formulae of quenched Chiral 
Perturbation Theory, which led to F n « 108 MeV. On the other hand, a method 
using only the zero mode contributions to the pseudoscalar correlators gave a 
lower value, close to the phenomenological F n . We conclude that these methods 
have the potential to determine Low Energy Constants of the chiral Lagrangian 
directly from the first principles of QCD, but precise values have to await the 
feasibility of simulations with dynamical, chiral quarks. 

As a step into this direction, we simulated dynamical overlap hypercube 
fermions in the Schwinger model using a Hybrid Monte Carlo algorithm with 
a simplified force, where we tested algorithmic requirements like reversibility [2]. 
The measurements confirm an extremely high level of locality. We also evaluated 
the chiral condensate based on the Dirac spectrum in the e-regime, and we found 
good agreement with low energy predictions at two very light fermion masses. 

References 

[1] W. Bietenholz and S. Shcheredin, Nucl. Phys. B 754 (2006) 17. 
[2] J. Volkholz, W. Bietenholz and S. Shcheredin, |hep-lat/0609003| 



35 



<t(600) as a Tetraquark Mesonium and the Pattern of 

Scalar Mesons 

Keh-Fei Liu 

Dept. of Physics and Astronomy, University of Kentucky, 
Lexington, KY 40506, USA 

Scalar mesons are not as well known as the pseudoscalar and vector mesons 
in terms of their SU(3) classification, their particle content and their spectrum. 
It has been known that there are too many experimental states for the qq nonet, 
possibly by a factor of two. We shall use lattice QCD calculations and combine 
with phenomenology to help understand the nature of the scalar mesons, their 
particle content and their SU(3) classification. 

A recent lattice calculation [T] with the overlap fermion on 16 3 x 28 and 
12 3 x 28 quenched lattices with the Iwasski gauge action at a = 0.2 fm has 
shown that the isovector scalar mesons with the ipip interpolation field has an 
unusual behavior as it approaches the chiral limit. Below the strange quark 
mass, its mass becomes flat and is almost independent of the quark mass, in 
contrast to other hadrons. This calculation which reaches a pion mass as low 
as 180 MeV is consistent with earlier findings with heavier quark masses [T]. 
Its chiral limit extrapolation indicates that ao(1450) is the qq meson and the 
ao(980) below it is not seen with the ipip interpolation field. It explains the 
experimental fact that Xq(1430) is basically degenerate with ao(1450). From the 
near degeneracy of a (1450), Kq(1430), and / (1500), we build a model for the 
mixing of / (1370), / (1500), and / (1710) through the annihilation of uu, (id and 
ss with a slight SU(3) breaking in their mixing and decays [3J. It is found that 
ratios of pseudoscalar meson decays can be well descried by this simple model. 
The interesting thing is that we found /o(1500) is almost a pure octet, / (1710) 
an almost pure glueball, and /o(1370) is mainly a singlet with ~ 10% glueball 
content. 

We also calculated tetraquark mesonium with 1=0 tttt interpolation field and 
found a state around 550 MeV in addition to the nir scattering states in the range 
where pion mass is lower than 250 MeV. Using the volume study of their spectral 
weight [3],[1], we verified that the lowest state around 300 MeV is the interacting 
7r7r state and the state at ~ 550 MeV is a one-particle state which we believe is 
the cr(600). Our present lattice calculation only confirms its existence. To get 
its precise mass and width, one needs to vary the lattice size to bring down the 
scattering state with one unit of lattice momentum to mix with it and observe 
how far the levels avoid each other. 

Based on the above two lattice calculations, we speculate on the following pat- 
tern for scalar mesons: ao(980), /o(980) and ft (800) form a tetraquark mesonium 
octet with the <r(600) as the singlet. On the other hand, ao(1450), 1^(1430), and 
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/o (1500) form a qq octet with / (1370) as its singlet and a almost pure glueball 
which is /o(1710). This needs additional experimental and lattice confirmation. 

This work is partially supported by the US DOE grant DE-FG05-84ER40154. 
I thank the organizers of the workshop on Lattice QCD, Chiral Perturbation 
Theory, and Hadron Phenomenology for giving me a chance to present this work. 
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We apply the formalism of covariant Baryon Chiral Perturbation Theory 
(BChPT) to the analysis of the moments of the Generalized Parton Distributions 
(GPDs) of the nucleon. We concentrate on the isovector n — 1 case and perfom a 
0{p A ) calculation to obtain the Generalized Form Factors A^^q 2 ), B^^q 2 ) and 
C% Q (q 2 ). These results can be used for chiral extrapolations of lattice data for 
GPDs. Of particular interest is the forward limit, where A2 Q (q 2 = 0) reduces to 
the averaged momentum fraction (x) u -d, a quantity known from phenomenology. 
Lattice data are also available for this moment. The figure shows two different 
chiral extrapolation functions. The dashed one indicates the influence of the well 
known leading chiral logarithm pQ, unable to reach the plateau observed in re- 
cent lattice simulations [2]. The full curve shows our 0(p 4 ) covariant BChPT 
result in the modified Infrared Regularization scheme [3]. We conclude [I] that 
a whole tower of quark mass dependent terms is required in order to understand 
(a;) u _rf(m 7r ), in contrast to other recent claims [5]. 
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ChPT for Generalized Parton Distributions (GPDs) 
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Generalized parton distributions provide a unified parametrization of hadron 
structure. They allow one to combine in the most efficent way the information 
obtained from various inclusive and exclusive reactions and permit to extract 
informations on the internal hadron structure which cannot be obtained directly 
by experiment. However, the extraction of GPDs from experiment is very non- 
trivial, especially to NLO and NNLO in a s , which is needed in view of the typical 
lowish Q 2 . Therefore, lattice calculations of moments of GPDs on the lattice are 
especially important for GPD physics pQ. The extrapolation of such results to 
physical quark masses requires a ChPT analysis of these moments. 
We performed such an analysis, both for the pion GPDs [2] and for the isosinglet 
nucleon GPDs [3]. We proceed by first constructing the leading order ChPT 
equivalents of the operators apearing in the OPE analysis of moments of GPDs. 
We then construct the next-to-leading contributions and calculate the 1-loop 
corrections to the leading ones. 

The results allow one to identify the leading analytic terms which apear and 
permit to relate some LECs to known physics. 
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The QCD vacuum as seen by overlap fermions 
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With the advent of overlap fermions, which are the cleanest known theoreti- 
cal description of lattice fermions implementing chiral symmetry, it has recently 
become possible to analyze the vacuum structure of Yang-Mills theory and the 
nature of the chiral phase transition of QCD from first principles using a fermionic 
approach. 

We report results based on (9(150) low lying overlap eigenmodes on quenched 
Liischer-Weisz configurations on various T=0 lattices that were obtained within 
the QCDSF-collaboration pQ. The analysis is also extended to dynamical fi- 
nite temperature configurations generated by the DIK-collaboration [2] using the 
clover- improved Wilson action. 

The distribution of the lowest modes in selected topological sectors and the 
spectral density is in good agreement with the finite- volume prediction of quench- 
ed chiral perturbation theory and Random Matrix Theory. In the vicinity of the 
finite-T phase transition the emergence of a gap in the spectrum is observed. 

Truncating the spectral decomposition of both the overlap-based topological 
charge density operator and the gluonic field strength tensor allows for a removal 
of fluctuations at the scale of the cutoff. Comparing the UV-filtered operators, we 
observe a strong correlation between the peaks of the topological charge density, 
the regions of high (anti-)selfduality and the domains of strong local chirality 
of the lowest localized eigenmodes. In the chiral-symmetry restored phase, on 
the other hand, the signal of local chirality and (anti-)selfduality is strongly sup- 
pressed. In contrast to the singular low-dimensional structure of the untruncated 
topological charge density, the UV-filtering reveals a structure of the T=0 vacuum 
much more reminiscent of semiclassical-like instanton pictures. 
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Contemporary lattice QCD calculations have highlighted the need to understand 
the dependence of QCD observables on the quark mass, in order to provide reliable 
extrapolation formulas. In addition to this motivation, studying the microscopic 
dynamical mechanisms involved in the transition into the chiral regime of QCD 
may shed light on the structure of the non-perturbative quark-quark interaction. 

In this work we investigate the non-perturbative quark-gluon dynamics in- 
volved in such a transition in the context of the Interacting Instanton Liquid 
Model (IILM) pp. By computing the nucleon and pion masses for a wide range 
of quark masses, we show that the IILM reproduces the existing lattice data 
for pion masses in the range ~ 500 — 600 MeV. Fitting the nucleon masses, 
using covariant Baryon ;\PT at order 0(p 4 ), we obtain the chiral coefficients in 
good agreement phenomenology. In particular the parameters C\ and C3 match 
well the values extracted from the analysis of the existing irN scattering data. 
On the other hand, the same chiral extrapolation performed using the available 
MILC data instead of the IILM data leads to chiral parameters C\ and C3 which 
are inconsistent with the experimental data. Such an inconsistency is probably 
due to the fact that the available MILC data have not been extrapolated to the 
continuum limit. We also show that the IILM agrees well also with Chiral Per- 
turbation Theory in the small quark mass regime. For example, we found that 
the spectral density of the Dirac operator calculated in the IILM agrees with the 
behavior expected from chiral perturbation theory for two flavors. Hence, we can 
argue that the IILM can be used to study the transition into the chiral regime 
of QCD. We identify a characteristic quark energy scale, m* = 80 MeV, which 
determines the boundary of the near-zero mode zone of the quark propagator in 
the instanton background. We argue that this scale governs the transition into 
the chiral regime and we discuss its physical interpretation. 
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Discussion: Lattice meets Chiral Perturbation Theory 
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I summarise the main points of the discussion held at the workshop: 

What is the status of lattice determinations of LECs? 

Even "simple" LECs such as F n are difficult to extract, signified by the large 
variation of results for this quantity. The compatibility of results from different 
simulations must be better understood, and an attempt to do this using recent 
data was presented pp. 

It was noted that the precision with which lattice determinations of LECs are 
quoted by some groups, is far too high, given the quality of the lattice data. 

Have we made contact between ChPT and lattice QCD? 

Often one finds that authors regard the presence of chiral logs in lattice data as 
the criterion for contact between lattice QCD and ChPT, since the coefficient 
of the chiral log is universal. It was pointed out [2] that this is not necessarily 
true, since this coefficient acquires regularisation-dependent contributions if chiral 
symmetry is not respected. 

In another contribution [3] the two-loop analysis for the quark mass depen- 
dence of F n [1] was compared to lattice data by NPLQCD [5] with pion masses 
as low as 300 MeV. While the two-loop formula describes the lattice data rea- 
sonably well, the claimed precision for the lattice determination of F n seemed 
exaggerated. 

What do ChPT people need from the lattice? 

More lattice data at smaller pion masses (m, < 300 MeV) must become available 
for more reliable determinations of LECs. In mesonic ChPT, the 0(p 4 ) LEC 
Lq is not well constrained by phenomenology and should be a target for lattice 
calculations. Also, there are few results for LECs in the baryonic sector. Here, 
the lattice should try to find ways to determine the so-called class II LECs, or 
symmetry breakers, since class I LECs having a dynamical origin (i.e. those mul- 
tiplying terms with derivatives) are more easily determined from phenomenology 

m 

It was suggested that a compilation of the status of determining LECs be 
made available, possibly in the form of a web-page. Some information can also 
be found in [Bj. 

What do lattice people need from lattice people? 

A standardised way of quoting results in order to facilitate comparison of results 
from several groups [7j. In particular, one standard choice of lattice scale (e.g. 
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ro) should be employed in addition to any other scale that collaboration might 
favour otherwise. 
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